Percolation and orientational ordering in systems of magnetic nanorods 
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' Based on Monte Carlo (MC) computer simulations we study the structure formation of a system of magnetic nanorods. Our 
model particles consist of fused spheres with permanent magnetic dipole moments, as inspired by recent experiments. The result- 
ing system behaves significantly different from a system of hard (non-magnetic) rods or magnetic rods with a single longitudinal 
dipole. In particular, we observe for the magnetic nanorods a significant decrease of the percolation threshold (as compared to 
"non-magnetic rods) at low densities, and a stabilization of the high-density nematic phase. Moreover, the percolation threshold 
is tunable by an external magnetic field. 
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1 Introduction 

Magnetic nanoparticle&i play nowadays a key role in a num- 
'ber of technological contexts from storage media and design 
of new functional materials^"-, to medical applications-^. In 
many cases, the magnetic particles are suspended in a non- 
magnetic carrier liquid such as water or oil, yielding a col- 
loidal suspension often called "ferrofluid". The most promi- 
nent ferrofluids involve spherical particles with typical sizes 
'of about 10 nanometers. These particles can be considered as 
I .single-domain ferromagnets; thus they have permanent mag- 
"netic dipole moments. The resulting anisotropic and long- 
jange dipole-dipole interactions between the spheres play an 
important role for their cooperative behavior. Indeed, already 
in zero field (and small packing fractions) the energetic pref- 
erence of head-to-tail configurations can lead to the forma- 
,tion of long chains and percolating (i.e., system-spanning) net- 
■works^, while at larger packing fractions, various ordered 
■Structures are observed—. Application of an external (static) 
"magnetic field on suspensions of magnetic spheres yields the 
formation of aligned chains and bundle s '^''^ . As a "byprod- 
uct" one observes drastic changes of the material properties, 
particularly the shear viscosity ("magnetoviscous effect") and 
the thermal conductivity '^. Thus, magnetic suspensions are 
.prime examples of complex fluids, whose internal structure, 
' phase behavior and dynamic rheological properties can be ef- 
ficiently controlled by external parameters ^^i^^ . 

Within this area of research, much effort is currently 
devoted to the synthesis and characterization of magnetic 
particles with anisotropic shape, examples being magnetic 
rods^— and cubes^. One key issue, e.g., in the case of 
rods, is to stabilize the chains (against shear) and thus, to 
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enhance the magnetoviscous effect observed in "simple" fer- 
rofluids. Moreover, anisotropic magnetic particles enable per 
definition a broader variety of self-assembled structures and 
patterns. However, compared to the case of magnetic dipolar 
spheres, the collective behavior of anisotropic magnetic parti- 
cles is much less understood. 

In the present article, we investigate structure formation 
phenomena in suspensions of a special class of magnetic 
nanorods (MNR) by Monte Carlo (MC) computer simulations. 
Our model MNRs consist of dipolar spheres which are perma- 
nently linked ("fused") into a stiff chain with internal head- 
to-tail alignement of the dipole moments. This model is in- 
spired by recent experiments of Birringer et ai— who used 
a self-assembly (aerosol-synthesis) process to produce mag- 
netic rods composed of aligned ironoxide spheres. The mag- 
netic field created by one of these rods is a superposition of 
the dipolar fields of the individual spherical particles. This is 
in contrast to earlier models of rod-like particles with single 
longitudinal (or transversal) dipole moment. Indeed, models 
of that type have already been intensely studied more than a 
decade ago, both by theory and by simulations—"—. At that 
time, interest was mainly driven by the desire to understand 
the phase behavior of polar liquid crystals. However, as we 
will see in the present study, the structural behavior of our 
new model MNRs differs strongly from that of single-dipole 
rods. 

One main goal of our study is to explore the percolation 
behavior of the MNRs. This is an interesting issue not only 
in the area of magnetic fluids (recall that already simple mag- 
netic spheres display pronounced chain and network forma- 
tion^). Indeed, the question of percolation is intensively dis- 
cussed also for general rod-like, colloidal particles, including 
prominent examples such as carbon nanotubes— "— . Research 
in this directions is generally driven by the desire to produce 
lightweight and, at the same time, highly connected materials 
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with mechanical, thermal or electrical properties that are en- 
hanced relative to their counterparts in the corresponding non- 
connected systems '^i'^-^^" —. Thus, a general aim of these ef- 
forts is explore conditions under which the percolation thresh- 
old, i.e., the volume fraction required for the formation of 
systems-spanning clusters, is decreased—. It has already been 
shown that such as decrease can be realized by an increase of 
the aspect ratio— as well as by other (interaction-related) fac- 
tors such as depletion effects^i^^i^. In the present study we 
show that the superpositioned magnetic interactions between 
our MNRs provide yet another mechanism that strongly en- 
hances the tendency for percolation. 

Further topics of this article concern the stability of long- 
range orientational order, as well as the impact of an external 
magnetic field—. In particular, we show that the MNR inter- 
actions tend to stabilize the nematic phase in sufficiently long 
MNRs. However, contrary to what has been found in a recent 
analytical study of non-magnetic rods^, the nematic transi- 
tion does not surpress the percolation. Finally, we briefly re- 
port on the impact of an external magnetic field. It turns out 
that already small magnetic fields lead to a significant decrease 
of the percolation threshold. 

The rest of this paper is organized as follows. In section IZTI 
we begin by defining our MNR model, followed by the details 
of the MC simulations in section IZ21 In section [3] we present 
our numerical results, starting with the effect of the dipolar 
interaction on the cluster formation. We then proceed to a 
discussion of the percolation transition and the appearance of 
global orientational order in sections Il!2] and l33] respectively. 
In section [3741 we show our results in presence of an external 
magnetic field. Our conclusions are summarized in section 2] 

2 Model and Simulations 

2.1 Model 

Our model system consists of stiff nanorods made out of sev- 
eral identical hard spheres of diameter a with a point dipole at 
their center (see Fig. [T]left). The positions of the spheres are 
fixed with respect to the center of mass of the rod, and the ori- 
entations of the dipoles are always aligned with the symmetry 
axis of the rod. We will use I to denote the number of mag- 
netic spheres composing a nanorod, which is also its length 
in units of CJ, therefore / = 1 corresponds to the dipolar hard 
sphere (DHS) case. The interaction energy between two rods 
; and j is then the sum of the pair interactions between their 
respective interaction sites (magnetic spheres). That is 



/, 






(1) 



where Rjj is the vector joining the centers of rods / and j with 
orientations Clj and D.j. Each rod is composed of Ij interaction 



sites, with r" denoting the position of site a of particle /, and 
m" is its magnetic moment. The interaction between sites is 
given by 



w{r,mf,mP): 



= if r < (7 

(mf-m^-3(mf •r)(m^-r)j/r3 ifr>c7 

(2) 



with r = rf — Tj . The magnetic moment has the same strength 
for each sphere and is given by m = \mf \ ~ jXrVsphMspi,, where 
Hr is the relative magnetic permeability of the solvent, and 
\'sph and Msph are the volume of the particles and the mag- 
netization, respectively. These parameters can, in principle, 
be extracted from experiments. Here we rather consider the 
reduced magnetic dipole moments m = \J^m~ /a^, where 
j8 = {IcbT)^^ with kg and T being Boltzmann's constant and 
the temperature, respectively. Common experimental values 
for m (at room temperature) are of the order 1 < '« < 10 '^'^^ . 
We also consider the effect of an constant external field (B) 
parallel to the z axis of the simulation box. The interaction 
energy between the external field and a MNR is 
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The coupling strength (relative to ksT) then follows as niB : 
j8mB. This expression suggests to define 



B=^/^B. 



(4) 



Experimental values for the magnetic fields are of the order 
of 0. 1 Tesla (implying B ~ 5 at room temperature with a = 
IQnm). MNR ferrofluids have been found to be susceptible 
already to small fields of B < lOmT (see Ref.— ), implying 
B<0.5. 

In the present study we also compare the structure of the 
magnetic nanorods consisting of / spheres (our MNRs) to 
the simpler, and well studied^J^"^, model of prolate hard 
spheroids with length 2a, width 2b, and a single longitudinal 
point dipole located at the center ("single dipole model"). We 
set the width to 2b = a. Thus, the shape and volume of a 
spheroid is comparable to that of a MNR if 2a — la. In or- 
der to compare the dipolar coupling strengths within the two 
models, we require that the pair interaction energies in the par- 
allel side-by-side configuration depicted in Fig.[T]are the same. 
This is 
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Using the previous condition, and the expressions for the in- 
teraction energies 
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and 
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we find that the "equivalent" dipole moment nig, as defined 
above, is given by 
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Fig. 1 Parallel side-by-side configurations of two MNRs (left) and 
two single-dipole spheroids (right). 



2.2 Monte Carlo simulations 

We carried out MC simulations in the NVT ensemble with pe- 
riodic boundary conditions for monodisperse systems of rods 
with lengths / — 4 and Z = 10, at several packing fractions 77 
(77 = INvsph/L^, where L is the length of the simulation box). 
Dipolar long range interactions were taken into account by us- 
ing Ewald summations with a conducting boundary—. 

Simulations were made with A^ = 1200 rods and A^ = 480 
rods for I = 4 and / = 10 respectively. In order to compare 
with the DHS case, we also simulate systems with / = 1 and 
A^ = 1000 particles. The points in phase space studied were 
averaged for 5 x 10"* to 2 x 10^ MC cycles, depending on the 
strength of the interaction. This relatively small number of 
steps is sufficient for a first, exploratory study, which is the 
aim of this paper. To speed up the sampling of the phase 
space cluster moves and inversion moves (changing the sign 
of the orientation of the particle) were used in addition to sin- 
gle particle moves. Cluster algorithms that speed up the sam- 
pling of the phase space in Monte Carlo simulations are widely 
used—"—. Here we define the clusters using a simple proxim- 
ity criteria approach, in which two particles are considered to 
be "bonded" (and part of the same cluster) if the nearest dis- 
tance between their surfaces is less than some value 5 (see 
below). 



The same cluster definition (which is based on a geometric, 
rather than an energetic criterion) is used to study the aggre- 
gation properties of the system. One important quantity in this 
context is the degree of polymerization ^ '1^^ , defined as 
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where A^^ / is the number of rods that belong to a cluster com- 
posed of more than one rod and A^ is the total number of rods. 
Second, we consider the percolation probability n defined as 
the probability of finding at least one "infinite" cluster, that is, 
a cluster connected to its own periodic images. 

For both properties, <I> and H, the choice of the parameter 
5 is in general arbitrary. From a physical point of view 5 may 
be thought of as the "hopping" distance for an excitation to 
go from the surface of one particle to another—. Depending 
on the microscopic nature of this excitation, the hopping can 
be relevant for, e.g., the thermal or electrical conductivity of 
the material. Here we consider 5 essentially as an adjustable 
parameter Specifically, following previous studies on perco- 
lation of nanorods 2'^i33,34 ^g ^^g ^ value for 8 which is small 
with respect to the dimensions of the particles (5 = 0.1 a). In 
addition, to get a somewhat less arbitrary measure of percola- 
tion, we also calculate the geometrical "critical" distance 8c- 
The latter is defined as the averaged minimum value of 5 for 
which an infinite cluster appear s ^''^^ . From a physical point 
of view, 8c can be interpreted as an inverse measure for the 
conductivity of the system—. 

For studying the structure of the system we computed some 
of the coefficients h'""' (r) of the expansion of the pair corre- 
lation function in terms of rotational invariants—. Here we 
present results for the projection— 
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h^'^ir) = 
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where p* = lvsph<^^ /l is the reduced density, P2 is the Legen- 
dre polynomial of degree 2, and Rij and 0,y are the center-to- 
center distance and the angle between the orientations of rods 
; and j. 

Finally, to investigate the degree of global orientational or- 
der we compute the degree of parallel ordering ("polariza- 
tion") and the nematic order parameters. The "polarization" 
is defined as 



n = T7 
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with d denoting the unit eigenvector associated with the 
largest eigenvalue of the matrix 
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where the / denotes the particle and the indexes k and / denote 
the cartesian component of the orientation vector The nematic 
order parameter (G2) is defined as the largest eigenvalue of Q. 
When applying an external field to the systems we measure 
additionally the magnetization of the system, defined as 
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3 Results 

3.1 Cluster formation 

Our first goal is to investigate the type of clusters formed in our 
MNR system at low packing fractions 77, as well as the depen- 
dence of the clustering on the interaction strength (m). The ap- 
pearance of clusters in systems of MNRs is already suggested 
by the well-studied system of DHS (corresponding to the case 
/ = 1 in our model). Indeed, DHS particles are known to form 
long, head-to-tail chains at packing fractions rj < 1 and dipole 
moments m > 2^. 

Taking the DHS system as a reference, we here present re- 
sults from MC simulations of rods with lengths / = 1 , 4, and 1 
and dipole moments m = 0, 1 .5, and 2.4, at a packing fraction 
7] = 0.0524. This packing fraction is comparable with typical 
experimental values (77 ss 0.01 — 0.02) for MNR systems—. 

Corresponding MC results for the degree of polymerization 
[see eq.|9l are given in Table[T] At m = (non-magnetic rods), 
the value of <I> is generally small and there is no formation of 
chains or other larger structures, regardless of the value of /. 
This is expected since the packing fraction considered here is 
far below the percolation threshold of that system (see Section 
3.2). Atm~ 1.5, about 30-40 percent of the rods are associ- 
ated into clusters [see Table [T| which contain, however, typi- 
cally less than three particles. Some representative snapshots 
of aggregates at m = 1 .5 and / = 4, 10 are shown in Figs.|2](a) 
and (b). It is seen that the rods already tend to form head-to- 
tail arrangements, despite of the still rather small interaction 
strength. Finally, at m = 2.4 essentially all rods are associated 
into clusters, as seen from the fact that <I> is close to one for all 
values of / [see Table [T|. As expected, these clusters contain 
on the average more rods than those at m = 1.5. The result- 
ing structures are illustrated by the snapshots in Figs. |2] (c) 
and (d). These snapshots also reveal that there are essentially 
three types of configurations which seem to be preferred by 
strongly coupled MNRs. Type I, which is the dominant one at 
1 = 4- [see Fig. 2 (c)], is given by head-to-tail chains or rings, 
similar to what is found in conventional DHS systems (/ = 1). 
Type II, which is seen in the / = 10 system, and, to less extent, 
also in the / = 4 system, is a side-by-side configuration where 
the dipoles point in opposite (antiparallel) direction. Finally, 
Type III, which seems to be particularly important for long 



MNRs (/ = 10), [see Fig. |2](d)] consists of a configuration 
where the dipoles are oriented along the same direction (par- 
allel), and the rods are close to side-by-side, but somewhat 
shifted against one another in longitudinal direction. 

Table 1 Degree of polymerization alrj = 0.0524. 
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Fig. 2 Examples of non-percolating clusters at a volume fraction 
77 = 0.0524. (a) m = 1.5 and / = 4. (b) m = 1.5 and / = 10. (c) 
m = 2.4 and / = 4. (d) m = 2.4 and / = 10. 

To better understand these configurations we take a closer 
look at the energetic landscape of interacting MNRs. A first 
useful insight is given by Fig. [3] which shows the interaction 
energy of a pair of MNRs as a function of the "joint" angle 7 
(in units of n). The latter is defined such that 7=1 and 7=0 
correspond to the Type I and Type II configurations, respec- 
tively. We see from Fig. [3] that for DHS (/ = 1) there is just 
one energy minimum with respect to 7, which corresponds to 
Type I. For / = 4 and Z = 10, on the other hand, there is a 
second pronounced minimum at 7 = 0. Moreover, this second 
minimum is separated from the first one by a large energy bar- 
rier which becomes the larger as / is increased. This suggests 
that both. Type I and Type II configurations are quite stable 
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against (thermal) fluctuations, a picture which is indeed con- 
firmed in the actual MC simulations. 



3^ 




Fig. 3 Interaction energy between two rods (i and j) composed of / 
magnetic spheres with m = 2.4, in Type I configuration as a function 
of the "joint" angle y, in units of ;r. y = 1 is the straight head-to-tail 
configuration and y = is the antiparallel configuration. 

Further, especially for longer MNRs (/ = 10) it is interest- 
ing to compare the energies related to Type II and Type III 
configurations. This is done in FigH] where we plot the inter- 
action energies of two adjacent rods (see sketch in the upper 
part of the figure) with parallel rod axes and either antipar- 
allel (Type II) or parallel (Type III) orientation of the dipole 
moments as function of the "longitudinal" distance d (i.e., the 
displacement of the rod centers along the direction of the rod 
axes). Note that we have plotted only the minima of the inter- 
action energies with respect to d. For the antiparallel configu- 
rations, these minima occur at d'"'" ^na with n being an inte- 
ger, whereas for the parallel configurations, d"™ = (n + 1 /2) <7 
corresponding to shifted (Type III) arrangements. From the 
numerical values of the two interaction energies we can see 
that both type of configurations are stable (in the sense that 
the energies are negative) for a broad range of displacements. 
Moreover, for larger values of d the energies related to Type III 
becomes even more attractive than those related to Type II. Fi- 
nally, for the specific values of / == 10 and the specific lateral 
distances considered in Fig. |4] the configurations with lowest 
energy are an antiparallel one with d = Q (i.e., no displacement 
at all) and a parallel one with d = 2.5(7. 

The presence of several types of preferred structures in 
strongly coupled MNR systems is also reflected by the cor- 
relation function h^^^{r) plotted in Fig.|5] At / = 1, one ob- 
serves the typical peaks at multiples of a particle diameter, 
reflecting the (head-to-tail) alignment of neighboring particles 
in the chains. The same type of (chain-like) structure is pre- 
ferred by the system with / = 4, as indicated by the large peak 
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Fig. 4 Left: Interaction energy minima between two adjacent 
parallel rods of length / = 10 with antiparallel (Type II) and parallel 
(Type III) dipole moments {m = 2.4). Right; Illustration of the 
configurations. 



of h^^'^{r) at r = 4c7. Note, however, that the / = 4 system also 
has a peak at r/a = 1 related to an antiparallel side-by-side 
(Type II) configuration. At Z = 10, the highest peaks occur at 
r = 1(7 and r w V^CT, whereas the (expected) peak at r = 10a 
is already to small to be distinguished from statistical noise. 
The large heights of the first two peaks reflect the presence of 
both. Type II and Type III configurations. 
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Fig. 5 The correlation function /i^^°(r) for MNRs with m = 2.4 at 
77 = 0.0524. 

To close this section, and to better acknowledge the spe- 
cific type of clusters formed in the present model of MNRs, 
we show in Fig. |6] two typical clusters formed in systems of 
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magnetic prolate spheroids with a single, longitudinal (point) 
dipole. For both aspect ratios considered in Fig. |6] the 
spheroids tend to align in an antiparallel side-by-side configu- 
rations, similar to the previously introduced Type II configura- 
tions with d=0. On the other hand, configurations of Type III 
are essentially non existent. The presence of antiparallel con- 
figurations in the single-dipole rods leads to quite compact 
clusters. This is consistent to findings in earlier MC studies^'* 
of such systems. 




Fig. 6 Sample clusters of hard spheroidal rods with a longitudinal 
point dipole at their center, with nie given by eq. ([8} with m = 2.4. 
The clusters are defined with a proximity criteria using 8 = 0.2b, 
and at a volume fraction r] = 0.0524. (a) a/b = 4. (b) a/b = 10. 



3.2 Percolation 

As a next step we want to understand the influence of the mag- 
netic interaction on the percolation properties of the MNRs. 
It is well known that percolation (i.e., formation of system- 
spanning clusters) already occurs in non-magnetic systems of 
prolate, hard-core particles, with the percolation density de- 
creasing upon increase of the aspect ratio of the particles ^^'^^ . 
In Fig. |7] we plot our present MC results for the percolation 
probability 11 of MNRs with I > 1 and various values of m as 
function of the packing fraction. All data in Fig. |7]have been 
obtained with a fixed value of the parameter 8 determining our 
cluster criterion (see section IZZt . 

Considering first the non-magnetic case, we see that both 
the systems with / = 4 and 10 do have a percolation tran- 
sition, as indicated by the steep increase of 11 from values 
close to zero to values close to one (upon increase of Tj). A 
precise determination of the percolation transition would re- 
quire calculations for different system sizes, which we did not 
perform in our present, more exploratory study. Neverthe- 
less, we can estimate from the data in Fig. [Tjthat the perco- 
lation thresholds at m = 0, using 5 = 0.1(7, are at 77/.='* « 0.21 
and Tj/.^'" « 0.14. To check the effect of 5 we performed 
calculations using 5 = 0.2(7 and obtained 7]/='* « 0.15 and 
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spherocylinders is important to remember two things: i) The 



Fig. 7 Percolation probability as a function of the volume fraction 
using d = O.Ict. 



length / of our rods is related to the aspect ratio of the sphe- 
rocylinderts via / = L/D + 1, where L is the length of the 
cylinder and D the diameter of the spherical cap of the sphe- 
rocylinder ii) At the same aspect ratio, the volume fraction 
of the multiple-sphere rod is lower than that of the sphero- 
cylinder. This means we have to multiply the volume fraction 
obtained for the multiple-sphere rod by a factor Vspcyi/vmsph^ 
where Vspcyi is the volume of the spherocylinder and v„„p/, is 
the volume of the multiple-sphere rod with / = L/D + 1 . Tak- 
ing into account these considerations, our results for m ~Q 
are in agreement with those reported for hard spherocylin- 
ders—. Coming back to the behaviour of MNRs, we note that 
the trend of longer rods percolating at lower densities, persists 
when we "switch on" the magnetic interactions, as revealed 
by the data for m =1.5. More importantly, we also see that 
the corresponding thresholds are smaller than those for m — 0, 
indicating that the magnetic interactions promote the perco- 
lation transition. This effect becomes even more pronounced 
for m = 2.4. For the latter case, we also see that the perco- 
lation thresholds seem to saturate for different rod length at 
TJ, w 0.03. 

Having obtained the percolation transition(s) for a fixed 
value of the cluster parameter 5, it is interesting to look at 
the percolation phenomenon from a less "biased" perspective, 
that is, without making any assumptions on particle distances 
within clusters. To this end we now consider the parameter 5,, 
defined in section 12.21 as the minimum value of 5 for which 
a percolating cluster appears in the system. Results for 5c as 
function of 77 are plotted in Fig. [H] As expected, all systems 
exhibit a decrease of 5, as the density increases. More im- 
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portantly, we see that already for m ~Q, the values of 5^ for 
the shorter rods (/ = 4) are consistently larger than those for 
/ = 10, indicating that the tendency for percolation is enhanced 
upon increasing the rod length. The same trend appears when 
we "switch on" the magnetic interactions, as shown by the 
data for m = 1.5 in Fig. [8j where the curves for fixed / are 
shifted to the left (i.e., to lower packing fractions) relative to 
the corresponding ones at m = 0. Thus, the magnetic interac- 
tions have the same supportive effect on the percolation. These 
findings are qualitatively consistent to the ones obtained for 
fixed 8 plotted in Fig.|2l 



to 




0.25 



Fig. 8 Geometrical "critical" distance (in units of a) vs. volume 
fraction for short (/ = 4) and long (/ = 10) rods with (m = 1.5) and 
without (m = 0) magnetic interaction. 

Finally, it is interesting to briefly compare the percolation 
behavior of MNRs system to that of systems of single-dipole 
spheroids. To this end, we consider the parameter 

A5r = — ^ X 100, (14) 

which measures, for fixed values of r\ and m, the value of 5c 
obtained for magnetic (MNR or single-dipole) particles rela- 
tive to the corresponding ones for m — {5'c). Results for 
A8r obtained for m = 2.4 and a low density are given in Ta- 
ble 12] In this context, the simulations have been started from 
an initial configuration with randomly positioned and oriented 
rods. The systems have then been equilibrated for ^^ 3 x 10^ 
MC steps. We see from Table |2] that for the MNR model 
the magnetic interaction produces negative values of A5r and 
thus, a decrease of 5c relative to the non-magnetic case, con- 
sistent with the results plotted in Fig. [8] For the single-dipole 
model, one the other hand, A5r is positive, meaning that the 
magnetic interactions rather hinder the percolation. 



Table 2 Comparison of A5r for different rod models with m - 
77 = 0.0524 and lUe given by eq. l[8j. 



:2.4, 



/ = a/b 



A5r(MNR) 



A5x(sing-dip) 



4 
10 



-87.8% 
-83.4% 



16.3% 
24.5% 



This remarkable difference in the behaviors of the two mod- 
els can be attributed to the different structure of the corre- 
sponding clusters. Whereas the MNRs tend to form open, 
elongated structures [see Fig. ID, the single-dipole spheroids 
rather tend to compact clusters with antiparallel local ordering 
[seeFig.|6l. 

3.3 Nematic ordering 

So far we have focused on the clustering and percolation be- 
havior of the magnetic nanorods. However, given the non- 
spherical, elongated shape of the MNRs it is tempting to con- 
sider also the appearance of long-ranged orientational (i.e., 
liquid-crystalline) ordering in our model. Indeed, it is well 
established that fluids of non-spherical particles can display 
entropy-driven phase transitions into nematic, smectic, and 
plastic phases, as the packing fraction is increased—. In par- 
ticular, the phase diagram of hard spherocylinders has been 
mapped out in detail by Bolhuis and Frenkel— who performed 
MC simulations for a broad range of length-to-breadth ratios 
L/D (where we recall that L/D = / — 1). According to their 
results, rods with L/D > 3.7 display both, a nematic and a 
smectic-A phase, the latter being characterized by the for- 
mation of layers in planes transversal to the global director 
Moreover, there are MC results for dipolar hard spherocylin- 
ders with L/D ~ 5 and a single longitudinal point dipolei^^. 
According to this study, the (longitudinal) dipole tends to 
destabilize the nematic phase. Instead, one observes an en- 
hanced tendency for smectic ordering. The stabilization of 
smectic-like ordering in these systems can be understood as a 
consequence of the strong energetic preference of antiparallel 
side-by-side configurations. Similar observations have been 
made for other model systems involving elongated particles 
with longitudinal dipole moment s ^^'^^ . 

In the present study we investigate the occurrence of liquid- 
crystalline (specifically, nematic) ordering via the orienta- 
tional order parameter Go defined in section 12.21 Results for 
G2 as function of the packing fraction T] are plotted in Fig.|9] 

For the case / = 4, the order parameter has negligible val- 
ues regardless of the dipole moment, indicating that there is 
no orientational ordering within the density range considered 
(77 < 0.25). Indeed, for the particular case m = 0, this result 
is consistent with the MC study of Bolhuis and Frenkel— who 
showed that hard spherocylinders with L/D + I = 4 do not 
present an isotropic to nematic phase transition. 
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Fig. 9 Nematic order parameter as a function of the volume fraction 
for several lengths and dipole moments. 



On the other hand, we find from Fig. |9] that the systems 
with / = 10 clearly do exhibit ordered phases, as reflected by 
the pronounced increase of G2 (upon increasing 77) towards 
values close to 1 . At the same time, the corresponding values 
of the "ferromagnetic" order parameter Gi defined in eq. ( fTTT ) 
are negligible, indicating that we indeed observe an ordering 
of the rod axes (rather than ordering of their dipole directions). 
Due to the relatively small systems sizes in our simulations, 
the G2-curves in Fig.|9]are rather rounded (in fact, a thorough 
determination of the packing fraction related to the ordering 
transition would require a systematic finite-size study). Nev- 
ertheless, we can extract two important findings from these 
curves. First, considering our model rods with / = 10 and 
m = 0, the behavior of the corresponding G2 appears to be 
consistent with the fact that the isotropic-nematic transition 
for pure hard spherocylinders (of aspect ratio L/D +1 = 10) 
occurs at 77 w 0.25-^ (indeed, taking into account the same 
considerations made in section[32]concerning the comparison 
of spherocylinders and multi-sphere rods, the packing fraction 
related to the onset of nematic ordering follows from Ref.— 
as 77 « 0.17). The second, and in the present context more 
important, finding concerns the influence of magnetic inter- 
actions. Indeed, we see that the functions 02(17) are shifted 
towards lower packing fractions when m is increased from 
zero. In other words, the nematic phase is stabilized (relative 
to the isotropic phase) due to magnetic interactions. This is 
is strong contrast to the behavior found in systems of dipolar 
spherocylinder s ^^-^^ (and other elongated particles with sin- 
gle, longitudinal dipoles— ), where the nematic state is rather 
suppressed. 

As an attempt to understand these differences, it is use- 



ful to consider the actual structure of the present MNR sys- 
tems within their nematic state. A simulation snapshot is 
shown in Fig.[TO](left). Closer inspection reveals that there are 
both, antiparallel (Type 11) configurations and shifted parallel 
(Type 111) configurations, consistent with our considerations 
concerning the energy landscape in section |3T| The "coexis- 
tence" of these configurations is also illustrated by the sketch 
in the right part of Fig. [TO] We suspect that it is particularly the 
presence of the shifted-parallel configurations, which could 
stabilize the nematic against the smectic-A phase formed in 
dense systems of (antiparallel oriented) single-dipolar elon- 
gated particles 



, 23.24 




Fig. 10 Left: Sample configuration of a system with / = 10, 
m = 2.4 and t] = 0. 1676 (G2 = 0.699). Right: Illustration of one 
way in which the magnetic interaction could stabilize the nematic 
phase. 



3.4 Effect of an external field 

In this last section we briefly discuss the impact of a homo- 
geneous, external magnetic field B = Bz on the structure of 
our MNR systems. We start by considering the norm of the 
field-induced magnetization (M=|M|, see eq. (fT3Tl). Results 
for the functions M(B) at two packing fractions, various rod 
lengths and a fixed reduced dipole moment of m = 1 .5 are 
plotted in Fig.[TT] For both values of 77 considered, the longer 
rods are seen to be more susceptible than the shorter rods, in 
the sense that the magnetization rises more sharply with the 
field and reaches earlier its saturation value. We also see that 
a decrease of B to zero leads to a vanishing of the magnetiza- 
tion in all cases considered, consistent with our result in sec- 
tion [53] that the Gi order parameter is zero (no spontaneous 
magnetization). 

A further interesting question is to which extent the external 
field influences the percolation phase transition. Some repre- 
sentative results in this context are shown in Fig. [12] Specifi- 
cally, in Fig. [T2l' a) we have plotted the percolation probability 
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Fig. 11 Magnetization of the feiTofluid with m = 1.5 as a function 
of the external field strength (6) for two volume fractions: (a) 
7] = 0.0524 and (b) rj = 0. 1676. 



Fig. 12 (a) Percolation probability vs. volume fraction for / = 10 
and m= 1.5 using 5 = O.lc. (b) Percolation probability vs. field 
strength for different values of I and 5. (c) The parameter Sc vs. 
field strength. In parts (b) and (c), the packing fraction rj = 0.0524. 



in a field of strength B = 5 as a function of 7] for rods of length 
/ = 10 and dipole moment m =1.5. Comparing this curve 
with the corresponding zero-field result [B = 0, also included 
in Fig. [T2l" a)l we find that the percolation threshold is shifted 
towards lower packing fractions. In other words, the mag- 
netic field supports the percolation transition in this system. 
In Fig. [T2l b) we present data for the percolation probability 
n as function of B at a fixed packing fraction. An increase 
of B enhances the percolation probability, consistent with our 
finding from Fig. [T2ja). One also sees from Fig. [T2l b) that 
the actual value of n in presence of a field strongly depends 
on the rod length (consistent with the zero-field situation) as 
well as on the the parameter 5 (determining our cluster crite- 
rion). To complement the discussion, we present in Fig. [T2l c) 
results for the "critical" distance 5c as function of the field 
and different rod length, /. For all lengths studied the value 
of 5c decreases upon an increase of B. This suggests that the 
percolation threshold is indeed reduced by the presence of an 
external field, regardless of the specific choice of 5. 

Finally, it is worth to briefly comment on the particle struc- 
tures observed in an external field. One general observation 
was that the structures of Type II (antiparallel local ordering) 
disappear, as expected from our results for the magnetization 
[see Fig. 11]. At the same time, structures of Type III become 
prevalent, along with a certain amount of Type I structures. 
The clusters that percolate under the effect of the external field 
do so in the direction of the field, that is, the system essen- 
tially percolates along one direction. This is in contrast to the 



zero-field situation, where the percolating clusters tend to span 
the simulation box in all three spatial directions. As a conse- 
quence, percolation in field leads to a strongly anisotropic net- 
work. This feature might become important e.g. for the design 
of systems with conductance anisotropy. 

4 Conclusions 

Based on MC computer simulations we have explored struc- 
ture formation phenomena in systems of model magnetic 
nanorods composed of / fused dipolar spheres. The parti- 
cle model has been inspired by recent experiments—, where 
MNRs have been produced via a self-assembly process. Com- 
pared to the experimental systems, we note that the longest 
rods considered in our simulations (/ = 10) have about half 
the length of their experimental counterparts (Z sa 24). This 
is because simulating longer rods would have required larger 
system sizes and thus much longer simulation (equilibration) 
times. Nevertheless, we are already in the realistic range. 
The same holds true for the range of reduced dipole mo- 
ments (m < 2.4) used in our simulations, which approaches 
the experimental value of m w 2.6 (for rods composed of 
iron nanospheres at room temperature-l^). Concerning the 
packing fractions, on the other hand, only the lowest value 
considered here (77 « 0.05) is close to current experiments^ 
(rj Ri 0.01 — 0.02). However, it is just the advantage of com- 
puter simulations that one can explore much larger parameter 
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ranges. 

Our simulation results generally reveal that the MNRs be- 
have strongly different not only from systems of individual 
magnetic spheres, but also from non-magnetic rods or rod-like 
particles with single dipoles. Most of our results refer to sys- 
tems in zero field. Compared to the case of magnetic spheres 
with central point dipoles (/ = !), we have found that the per- 
colation threshold for / > 1 is lowered towards significantly 
lower packing fractions. At the same time, the percolation 
thresholds are also much lower than those of non-magnetic 
rods of comparable lengths. These findings indicates that 
the MNRs could be promising candidates as building blocks 
of lightweight nanocomposites, i.e., connected materials with 
novel mechanical and conductivity properties. Moreover, for 
sufficiently large / and densities beyond the percolation thresh- 
old, the MNRs display an isotropic-to-nematic transition, but 
no long-range ferromagnetic ordering, which is again in con- 
trast to magnetic spheres. These differences regarding the per- 
colation and ordering behavior can be explained, on a micro- 
scopic level, from the fact that the present MNRs display a 
larger variety of cluster structures. In particular, in addition to 
the usual head-to-tail ordering of dipolar spheres, we observe 
formation of local antiparallel ordering, as well as of parallel 
side-by-side ordering with shifted positions. 

The presence of the latter type of clusters also represents 
one main difference to systems of rod-like particles with single 
(central, longitudinal) dipole moments. Indeed, this simpler 
model mainly exhibits compact clusters with local antiparallel 
ordering—. As a consequence, these single-dipolar rods rather 
exhibit smectic (instead of nematic) phasesSJ^^. From a more 
general perspective, our simulations results therefore suggest 
that the distribution of dipoles within the nanorods is of cru- 
cial importance for both, their percolation behavior and for the 
nature of ordered phases. 

Finally, we have briefly discussed the impact of a static 
magnetic field. Our main conclusion in this context is that 
already relatively weak fields (with a strength of about 2 mT) 
can again significantly lower the percolation threshold com- 
pared to the zero-field situation. 

In the present paper we were concerned with the 
(interaction- or field-driven) structural behavior of the MNRs 
in thermal equilibrium. From an experimental point of view, 
it would be very interesting to explore the impact of these 
structural phenomena on the single-particle dynamics (e.g., 
the translational and rotational diffusion) and the rheological 
properties of the suspension (in particular its shear viscosity). 
Clearly, an investigation of these phenomena requires to sup- 
plement the MC simulations employed in the present work by 
computer simulations targeting the time-dependent behavior, 
such as Brownian dynamics. Work in these directions is in 
progress. 
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